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Abstract 

In this paper we investigate the fermionic condensate and the renormalized vacuum ex- 
pectation value (VEV) of the energy-momentum tensor for a massive fermionic field induced 
by a flat boundary in the cosmic string spacetime. In this analysis we assume that the 
field operator obeys MIT bag boundary condition on the boundary. We explicitly decom- 
pose the VEVs into the boundary-free and boundary-induced parts. General formulas are 
provided for both parts which are valid for any value of the parameter associated with the 
cosmic string. For a massless field, the boundary-free part in the fermionic condensate and 
the boundary-induced part in the energy-momentum tensor vanish. For a massive field the 
radial stress is equal to the energy density for both boundary-free and boundary-induced 
parts. The boundary-induced part in the stress along the axis of the cosmic string vanishes. 
The total energy density is negative everywhere, whereas the effective pressure along the 
azimuthal direction is positive near the boundary and negative near the cosmic string. We 
show that for points away from the boundary, the boundary-induced parts in the fermionic 
condensate and in the VEV of the energy-momentum tensor vanish on the string. 



PACS numbers: 03.70.+k, 98.80.Cq, 11.27.+d 



1 Introduction 

Non-trivial structure of the vacuum state is among the most profound predictions of quantum 
field theory. The properties of the vacuum are manifested in its response to different external 
perturbations. The imposition of boundary conditions on a quantum field is among the simplest 
models for the perturbation. The boundary conditions may be induced by real material media 
or represent non-trivial topology of the space. In both cases they modify the structure of the 
quantum vacuum thus affecting the vacuum expectation values (VEVs) of physical observables. 
This phenomenon is known as the Casimir effect (for reviews see Ref. [I]). In the present paper 
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we consider combined effects from boundaries and non-trivial topology on the properties of the 
fermionic vacuum in the cosmic string spacetime. 

The cosmic strings are predicted to be formed within the framework of grand unified theo- 
ries as a result of symmetry braking phase transitions [2j. This type of phase transitions have 
several cosmological consequences and provide an important link between particle physics and 
cosmology. Although the recent observational data on the cosmic microwave background radi- 
ation have ruled out cosmic strings as the primary source for primordial density perturbations, 
they are still candidates for the generation of a number of interesting physical effects such as the 
generation of gravitational waves, gamma ray bursts, high energy cosmic rays, and enhanced 
baryon number violation. Recently, cosmic strings attracted a renewed interest partly because 
a variant of their formation mechanism is proposed in the framework of brane inflation [3] (for 
a review of the existence of cosmic strings in string theory models see [1]). 

At distances larger than the core radius, the spacetime geometry for a straight cosmic string 
is well approximated by a flat spacetime with a planar angle deficit. In quantum field theory, 
the corresponding non-trivial topology gives rise to the polarization of the quantum vacuum. 
This problem has been studied in many papers for scalar, fermionic and electromagnetic fields 
[5]- [31]. The vacuum polarization by a cosmic string in various curved backgrounds has been 
discussed in [32j [33j [34J, [35]. Another type of vacuum polarization arises by the presence of 
boundaries on which the field operator obeys some prescribed boundary condition. The analysis 
of the boundary-induced quantum effects in the cosmic string spacetime have been developed 
for scalar [36l [37] , fermionic [381 EH] > an d electromagnetic fields [37J [301 [31] for the geometry of 
a coaxial cylindrical boundary. From the point of view of the physics in the exterior region, the 
latter can be considered as a simplified model for the nontrivial core of the cosmic string. The 
Casimir force for massless scalar fields subject to Dirichlet and Neumann boundary conditions 
in the setting of the conical piston has been recently investigated in [32] • The Casimir densities 
for scalar and electromagnetic fields induced by flat boundaries perpendicular to the string have 
been considered in [431 144j . The topological effects in the generalized cosmic string geometry 
with a compact dimension along its axis are investigated in [35]. Combined effects of topology 
and boundaries on the Casimir-Polder force acting on a polarizable microparticle situated near 
the cosmic string have been studied in [36]. In particular, it has been shown that due to the 
nontrivial topology this force can be repulsive. 

Continuing in this line of investigations, in the present paper we consider the polarization 
of the fermionic vacuum by a flat boundary perpendicular to the string axis on which the field 
operator is constrained by MIT bag boundary condition. We evaluate the fermionic condensate 
(FC) and the VEV of the energy-momentum tensor. These quantities are among the most 
important characteristics of the vacuum state. Although the corresponding operators are local, 
they carry an important information about the global properties of the background spacetime. 
In addition to describing the physical structure of the quantum field at a given point, the 
energy-momentum tensor acts as the source in Einstein's equations and, therefore, it plays an 
important role in modelling a self-consistent dynamics which involves the gravitational field. 
The FC is crucial in the models of dynamical chiral symmetry breaking and in considerations of 
the stability of fermionic vacuum. As it will be shown below, in the problem under consideration 
all calculations can be performed in a closed form and it constitutes an example in which the 
topological and boundary-induced polarizations of the vacuum can be separated in different 
contributions. Note that a conical space appears as an effective background geometry in the 
long- wavelength description of certain condensed matter systems such as crystals, liquid crystals 
and quantum liquids (37J. The results described in this paper may shed light on the features of 
the boundary-induced effects in this type of medium with a conical defect. 

The rest of the paper is organized as follows. In the next section we describe the geometry 
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of the problem and a complete set of positive- and negative-energy four-spinors obeying the bag 
boundary condition on the plate is constructed. In section [3] these mode functions are applied 
for the evaluation of the FC. Similar calculations for the VEV of the energy-momentum tensor 
are presented in section [H The main results of the paper are summarized in section \E\ The 
technical details for the evaluation of the energy density and the azimuthal stress are described 
in Appendix lAl 



2 Geometry of the problem and mode functions 

We consider the background spacetime corresponding to an infinitely long straight cosmic string. 
Considering the string along the z-axis, by using cylindrical coordinates the line element asso- 
ciated with this spacetime reads: 

ds 2 = dt 2 - dr 2 - r 2 d<\> 2 - dz 2 , (2.1) 

where ^ 4> ^ 4>o = 2vr / q and the spatial points (r, (f>, z) and (r, <p + 4>o, z) are to be identified. 
The line element (|2. 1[) has been derived in [48j in the weak-field approximation. In this case, the 
planar angle deficit is related to the mass per unit length of the string, hq, by 2tt — 0o = 8irGfiQ, 
where G is the Newton gravitational constant. The validity of the line element (|2.ip has been 
extended beyond linear perturbation theory by several authors [39] (see also [2]). In this case, 
the parameter q need not to be close to 1. This is the case also in condensed matter systems 
with conical defects. For example, in graphitic cones one has q = 6/5,3/2,2,3,6. All this type 
of cones have been experimentally observed [50J. 

The dynamics of a massive fermionic field in curved spacetime with the metric tensor g^ v is 
governed by the Dirac equation 

i 7 ^V M V -m4 = 0, (2.2) 
with the covariant derivative operator 

V M = <9 M + r„ , r M = 7^7^/4 • (2.3) 

Here 7^ are the Dirac matrices in curved spacetime and is the spin connection. Here and 
below, \i = 0, 1, 2, 3 correspond to the coordinates t, r, (J), z, respectively, and the semicolon means 
the standard covariant derivative for vector fields. 

In the present paper we are interested in quantum vacuum effects induced by a flat boundary 
located at z = 0. We assume that on the boundary the field operator obeys MIT bag boundary 
condition: 

(l + ^n M )V = , z = 0, (2.4) 

where is the outward-pointing normal to the boundary. In the discussion below we will 
consider the region z > with = (0,0,0,-1). For the evaluation of the VEVs we will use 
the direct mode summation procedure. In this approach one needs to have the complete set of 
the eigenfunctions satisfying boundary condition (12.41) . 

In the discussion below the gamma matrices will be taken in the Dirac representation: 

with ft , I = 1, 2, 3, being the Pauli matrices in cylindrical coordinates with a planar angle deficit: 
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and /3 3 = diag(l,— 1). We have the standard anticommutation relation (3 l (3 n + f3 n {3 1 = —2g ln . 
For the spin connection one gets 

r„ = ^nWsl- (2.6) 

By taking into account that 7^1^ = (1 — q)^ 1 /(2r), the Dirac equation takes the form 



7 M <9^ + -=-^7V + »mV> = 0. (2.7) 
2r 



First we consider the positive-energy solutions of this equation obeying the boundary condition 
(El). 



2.1 Positive-energy modes 

For the positive-energy solutions the time dependence has the form e~ lujt . Decomposing the 
four-spinor ip into two-component spinors, 

,,s, 

from Eq. (|2,7p for the upper and lower components we find the set of equations 

pldl+ ^r pl ) v ~ i (w + m) x = °' 

tfdt+^-A X - i (w - m) y> = 0, (2.9) 



From (12.911 the following equation is obtained for the upper component: 



r 



4^2 



where i,j = 1, 2, 3. 

Further decomposing the two-spinor ip as 



9 %3 d i d j +-d 1 + - — -p v pd2- ™ A :' +uj 2 -m 2 



<p = 0, (2.10) 



we can present the solutions for the separate components in the form 

<Pl = Mr) (Cf ie lkz + Cf 2 e~ lkz ) e^'^, 1 = 1,2, 



(2.12) 



with ^ k < oo and ni = 0, ±1, ±2, .... Substituting (|2.12p into (|2.10p . we can see that for 
the radial functions the Bessel equation is obtained with the solutions regular at the string 
fi(r) = «/ft(Ar), where 

A = \qni - (-l) l (q - 1)/2|, (2.13) 
and A = \/uj 2 — k 2 — m 2 . By using (|2.12p . for the components of the lower two-spinor, 

X=( Xl ), (2-14) 



. X2 

we get 



Xi = JbAXt) [B«le ikz + B%- lkz ) e ^- u \ (2.15) 
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with the relations n 2 = Hi + 1, fa = fa + e m , and we have defined 



l, m ^ o 

-1, ni < 



(2.16) 



The coefficients in f|2. 15[) are defined by the expressions 



B 



B 



(1) 




^2 - " 




+ i\e ni 


v \ - 


oj + m 




oj + m 


(2) 




> n ip2 — 






V l ~ 


oj + m 




oj + m 



(2.17) 



Hence, the positive-energy solution to the Dirac equation has the form (n% = n) 



( 



(+) 



rvX) „ikz I z"f 0-) ik 

U ipl e + °^2 e 

dC 1 ) ifcz I ffWp-ik 

n <pl e + n ip2 e 



>(2) ife2 



>(2) n —ikz 
> ip2 



J^ +eB (Ar)e^ 



J,(qn4>—u>t) 



(2.18) 



where 



/3 = |g(n + l/2)-l/2|, 



(2.19) 



and o" is a collective notation for the set of quantum numbers specifying the solution (see below). 
We can see that the four-spinor is an eigenfunction of the projection of the total momentum 
along the cosmic string: 



with 



Note that €j = e r 



j = n + 1/2. 



(2.20) 
(2.21) 



Now we impose the boundary condition (|2.4ft on the four-spinor (I2.18|) . From this condition 
the following two equations are obtained for the upper and lower components: 



¥>i _ iXl =0, tf2 + iX2 = at z = 0. 



(2.22) 



By using these relations, the coefficients in (|2.18p are expressed in terms of and C^/ as 



C 



C 



(2) 
<pl 

(2) 
ip2 



JlL{[k 2 + m(oj + m)] C$ + {uj + m){m + ik) C$ } , 
- ^ | (w + m) ( m _ iA .) crW + [k 2 + m(uj + m)] C$ } , 



(2.23) 



and 



D 



B 



B 



B 



(l) 
<pi 

(i) 

(2) 
Vl 

(2) 
V 2 



m r (l) 
m 



m 



A: 



+ i C, 



,(i) 

>2 ' 



fc 

A L 
A L 



m 



,(i) 



(2.24) 
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One of the coefficients and is determined from the normalization condition obeyed 
by the wave function and, hence, one coefficient remains arbitrary. In order to determine the 
latter some additional condition should be imposed. A similar situation to find fermionic wave 
function was present for the problem with cylindrical boundary [38]. However, as it can be seen, 
the additional condition used there is not appropriate for the present case. 
In order to specify the second constant we impose the condition 



R (2) /R (2)_ r W /r W 
By taking into account Eq. (|2.24p . this leads to the following relation 



C 



(i) 



where we have defined 



ip2 

oj + sX 



%Kjc C, 



(1) 



±1. 



(2.25) 

(2.26) 
(2.27) 



k — im 

Note that we have the relation k s k*_ s = 1. 

With the relation (I2.26D all coefficients are expressed in terms of c5 \ Introducing the 
notations 



9±(z) 



e ikz ±iK,e~ ikz , 



(u ± m) f±(z) + s\f T (z), 
the positive-energy solutions are written in the form 



(2.28) 



/ f + (z)MXr)e-^/ 2 \ 



(+) 



c 



(1) 



Jq<P/2 



lg_(z)MXr)e-^/ 2 
V -ie r sf-(z)Jp +ej (Xr)e 1 ^ 2 J 



where j = ±1/2, ±3/2, . . ., and 

P = Px = Q\j\-e j /2. 
Now the set of quantum numbers is specified to a = (A, k, j, s) and 



(2.29) 

(2.30) 
(2.31) 



VX 2 + k 2 

It can be checked that the mode functions (|2.29[) are orthogonal: f d 3 x \/\g\ip^ + Vi' = 
for a a'. The remained coefficient C$ is found from the normalization condition 

(2.32) 



<l i x^/\g\^ + ^ ) = 5(X' - X)5(k' - k)5 fj 5 s , & 
By taking into account that 



drrJp{Xr)J p {X'r) = -5{X - A'), 
A 



we find 



IC 



(1)|2 



k 2 X 



(2.33) 



(2.34) 



^ ' 8TT<f) Uj(u + sX) ' 

It can be seen that, the positive-energy mode functions (j2.29p . with the normalization coefficient 
given by (|2.34p . are invariant, up to a phase, with respect to the transformation (k,s) —> 
(-k,-s). 
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2.2 Negative-energy modes 



The time-dependence of the negative-energy mode function is in the form e lujt . For the lower 
two-spinor, \t we obtain the same equation as in (|2.10p . For the corresponding upper and lower 
components one finds 



XI 



J Pl (\r)(C^}e ikz + C^e-^Je^^, I = 1,2, 



'A-2 



(2.35) 



with the notations (|2.13p . Substituting these expressions into (|2.9p . for the upper components 
we find the expressions 



m = {B% ikz + B%- ikz )J Pl (\r)e i{ - qniM \ 



(2.36) 



and the relations n<i = n\ + 1, fii = /3i + e ni . Denoting ni = n, the coefficients in (|2,36p are 
given by the expressions 



B 



B 



( ¥ = 




' J °X2 ~ 






w + m 


u + m 


(2) 


frr {2) ,-Af r (1) 


^X2 ~ 




xi - 


cj + m 


u + m 



(2.37) 



As a result, he negative energy modes are presented as 



(-) 



with p defined in (l2TT9j) . 



/" (B$e ikz + B { ^e-' ikz )J l3 (\r) \ 



i(qncj>+ut) 



(2.38) 



From the boundary conditions (|2.22p the coefficients are expressed in terms of eg and cty '■ 



C 



C 



(2) 
Xl 

(2) 
X2 



^g- { [k 2 + m (w + m)] + (w + m) (m + ifc) eg } , 

~5x {( w + m ) ( m - iA: ) cji + [ fc2 + m ( w + m )] c®} » 



(2.39) 



and 



I? 



(i) 



-°X2 



/m \ (i) m a) 



B 



B 



(2) 
Xl 

(2) 
X2 



En 
?A 



Co' 



+ (m + ik) C x2 
(m - ik) eg + wCg 



(i) 



(2.40) 



Similar to the case of positive-energy modes, an additional condition should be imposed. As 
such a condition we impose 

^X2/^X1 _ ~°X2/ U X1 

Taking into account (|2.4U|) . this gives 



(2)/ R (2) _ _ n Q-)inV) (2.41) 



C 



(i) 

X2 



2K s C.g, S 



±1, 



(2.42) 
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with k s defined in (|2.27|) . 

As a result, the negative-energy mode functions are written in the form 



/ 



(-) 



C 



(i) 



z)J /3 (Ar)e-^/ 2 \ 
ie iS /_(z)^ +ej (Ar)e^/ 2 
/ + (z)J^(Ar)e-^/ 2 
V %^-g + {z)J^{\r)e^ ) 



J(qj<t>+ut) 



(2.43) 



with the same notations as in (|2.29[) . The function is an eigenfunction for the projection of the 



total angular momentum on the z-axis, J^i- = qjipa ' , and we have introduced the related 
quantum number j = n + 1/2. It can be checked that the modes (|2.43j) obey the orthogonality 
relations J d?x^J\g\il) , ipi = for a ^ a' . From the normalization condition we can see that 

, a (D|2 



(") 



C$\ 2 . Note that we have a correspondence between the negative- and positive-energy 



modes given by the relation (up to a phase) 



-iuit 



iujt 



(2.44) 



where in cylindrical coordinates 7 5 = irj ^^/ 2 ^. This relation corresponds to the PCT trans- 
formation of four-spinors. 

The wave functions obtained in this section can be used for the investigation of various 
quantum electrodynamics effects around the cosmic string in the presence of boundary involving 
electrons and positrons (for example, bremsstrahlung and electron-positron pair production by 
a photon). In what follows we use these functions for the evaluation of the FC and the VEV of 
the energy- momentum tensor. 



3 Fermionic condensate 

Having the complete set of mode functions for the fermionic field, we can evaluate the effects of 
the boundary on the FC, (0|^^|0), where |0) is the vacuum state in the geometry of a cosmic 
string with flat boundary and ip = Tp +/ y° is the Dirac adjoint. Expanding the field operator 
in terms of the complete set {ipi- \tp^} and by using the definition of the vacuum state, the 
following formula is obtained for the condensate: 

(0|^|0)=^^-)^- ) , (3.1) 

a 

where 



oo roc 



£= E E / dX dk - ( 3 - 2 ) 

Substituting (|2.43[) into (|3.ip . the expression obtained is divergent and some regularization 
procedure is needed. We will assume that a cutoff function is present without explicitly writing 
it. The specific form of this function will not be important for the further discussion. 
By making use of the expression (|2.43j) for the mode functions and the relation 

£ j|(Ar)= £ J| +e ,(Ar)=^[j 2 ._ 1/2 (Ar) + J 2 +1/2 (Ar)], (3.3) 

j=±l/2,... j=±l/2,... j 

the FC is presented in the form 

(0|^|0) = (0|^V|0)s + W)b- (3.4) 



8 



Here the separate terms on the right-hand side are given by the expressions 



[* OO Z 1 OO \ 

«#)>. = / d\[ *-[jJ_ 1/a (Ar) + J5 +1/a (Ar)], (3.5) 



and 



■ /* OO OO \ 

W)b = ^E/ dA / ^-[^-i/ 2 (Ar) + 4 +1/2 (Ar)] 

im) e 2ifcz — (k + im) e~ 2ikz ] , (3.6) 



[(* 



where and in what follows we identify 



E = 

i j=l/2,3/2 



First we consider the part (0|V^|0) S . This part corresponds to the FC in the boundary-free 
cosmic string spacetime (sec |38j). The geometry of cosmic string is flat outside the string core 
and, hence, the renormalization is reduced to the substraction of the corresponding quantity in 
Minkowski spacetime: 

W)s,ren = (<#^|0) 8 - (0|^|0) M - (3.8) 



For the further transformation, we use the representation 

12 



After the substitution of (|3.9p into f)3.5j) and changing the order of integrations, the integral 
over A is expressed in terms of the modified Bessel function of the first kind, I v {z) (see |51|). 
Integrating over k and introducing a new integration variable y = r 2 /(2s 2 ), the boundary-free 
part is presented in the form 

Q m r .7_-«-m 2 r 2 / 



(0\W\0) s = -JL^J dye-y~ mr / 2 yi(q,y), (3.10) 
with the notation 

z(.q, v) = £ i^iMv) + i<u+i/2(v)] • (3-ii) 

3 

For the series (|3.1ip we use the formula 

v 

Z(q,y) = ~y"'(-l) l cos(Trl/q)e ycoB( - 27rl M 
H 1=0 

+ 2 - cos (W2) f d ^ff^^K -y^^ (3.12) 
7r Jo cosh(gx) — cos{qir) 

where p is an integer number defined by 2p ^ q < 2p + 2 and the prime on the sign of the 
sum means that the term I = should be taken with the coefficient 1/2. Formula (|3.12p is 
obtained as a special case of the more general formula derived in |39j. By taking into account 
that 1(1, y) = e y , we see that the contribution of the term / = in (|3,12p into {0\ipip\0) s coincides 
with the FC in Minkowski spacetime. The renormalized value of the FC in the boundary-free 
cosmic string spacetime is obtained subtracting this term. 
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Combining (|3.10|) and (|3.12|) . after the integration over y, for the renormalized FC in the 
boundary-free cosmic string spacetime one findd_| 



2mr 



2q 

H cos (qn/2) 

TT 



(- 1)' cos (tt/ /q) h {2mr sin(vr/ /g) ) 
sinh (qx) sinh(x)/i (2mr cosh x) 



l=i 

•oo 



dx- 



cosh(2qx) — cos(qir) 



where we have defined the function 



Mv) = y- v K v {y). 



(3.13) 



(3.14) 



Note that for this function one has the relation d y f u (y) = —yf u+ i(y). As it is seen from (|3.13|) . 
for a massless field the renormalized FC vanishes in the boundary-free cosmic string spacetime. 

Now we turn to the term (ipi^)b m (|3.4p which is induced by the boundary. For points away 
from the boundary, this term is finite and the cutoff function, implicitly assumed before, can be 
safely removed. For the further transformations it is convenient to write the boundary-induced 
term in the FC in the form 



l-OO [-CO 1 

W>b = A (2m + 3 Z ) V / d\\ dk- 

x cos(2kz)[J^_ 1/2 (Xr) + jJ +1/2 (Ar)]. 



(3.15) 



Similar to the boundary-free case, by using the representation (|3.9|) and changing the order 
of integrations, the integrals over k and A are explicitly evaluated. Changing the integration 
variable to y = r 2 /(2s 2 ), we get 



(2m + d z ) I dye 



-(2z 2 /r 2 +l)y-m 2 r 2 /( 



(3.16) 



with the notation (|3.1ip . As it is seen, the boundary- induced part in the FC does not vanish 
for a massless field. First we consider this case. 

For a massless field the integral in (|3.16p is evaluated by using the formula from [51] and 
after that the series over j is reduced to the sum of a geometric progression. In this way we find 



4:7r 2 rz 2 u 2q — 1 



1 + 



qz 



u 2q + l 



Vr 2 + z 2 u 2 i - 1 



(3.17) 



with the notation 



u = z/r + ^1 + z 2 /r 2 . (3.18) 

As it is seen, for a massless field the boundary induced part is always negative. By taking into 
account that in this case the boundary- free part vanishes, we conclude that the same is the case 
for the total FC. For q = 1, Eq. (|3.17p is reduced to the known result for the plate in Minkowski 
spacetime: 

' (3.19) 



W> b M) 



Air 2 z 3 ' 



(V^b^ i s t ne leading term in the asymptotic expansion for (ipip)b m the limit z — > 0. Assuming 
z <C r, we have 



1 



Air 2 z 3 



1 



(g -1)(V + 17 ) (z\± 
360 



(3.20) 



An alternative integral representation for (V'V^s.ren is given in [38] 
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For large values of z, z > r, the leading term has the form 



Note that for q > 1 the boundary-induced part vanishes on the string. 

Now we return to the case of a massive field. In order to provide an expression of the 
condensate more convenient for numerical calculations, we use the formula (I3.12p . After the 
substitution into (|3. 16|) . the integral over y is expressed in terms of the Macdonald function. As 
a result for the boundary-induced part in the FC we find 



2m 3 

7T 2 



(—1) cos {nl/q) F(2mz, 2mr sin(irl/q)) 

1=0 

2 f°° sinh (qx) sinhx 

H — g cos ((77172) / ax — — r -. — rF(2mz, 2mr coshx 

7T ,/ cosh(2gx) — cos(g7r) 



where 



F{x,y) = /i(Vx 2 + y 2 )-x/ 2 ( V / x 2 + 2 / 2 ). (3.23) 



In this formula, the I = term corresponds to the FC induced by a boundary in Minkowski 
spacetime: 

2 3 

(^)t M) = ^{2m + d z )h{2mz) = 2^F(2mz,0). (3.24) 
7r z ir z 

Near the boundary this term dominates in (|3.22p . The remained part corresponds to the con- 
tribution induced by the nontrivial topology of the string. For r / the latter is finite on 
the boundary, z = 0. In order to find the asymptotic expression of the boundary-induced 
FC for points near the string it is convenient to use the representation (|3.16p . The dominant 
contribution comes from the term with j = 1/2 in (|3.11|) and to the leading order one finds 

H 2(9- 1 )/ 2 7r 2 r((g + l)/2) V ' 

As it is seen, similar to the case of a massless field, for q > 1 and z ^ the boundary-induced 
part of the FC vanishes on the string. At large distances from the boundary, mz 3> 1, the 
boundary-induced part is exponentially suppressed. 

Having the separate parts, for the renormalized FC we get: 

(V^ren = W)s,ren + W)b- (3-26) 

In figure [U for the cosmic string with q = 3, we have plotted the quantity (ipip) ien /m 3 as a 
function of distances from the boundary and from the string in units of the Compton wavelength 
for the fermionic particle. Near the boundary, the boundary induced part dominates and the 
FC is negative. For points near the string, the FC is dominated by the boundary-free part and 
it is positive. 



4 Energy- momentum tensor 

In this section we consider another important characteristic of the fermionic vacuum, the VEV 
of the energy-momentum tensor. This VEV can be evaluated by using the mode sum formula 

(oivio) = \YM )( ^ ) ^ V ^ ) ~ (Vo^h*)^] , (4.1) 
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where the brackets enclosing the indices mean the simmetrization. Inserting the expression for 
the negative-energy mode functions, similar to the case of the FC, the VEV is presented in the 
decomposed form: 

(0\T ia/ \0) = (0\T flv \0) B + (T flv ) h , (4.2) 

where the first and second terms on the right-hand side correspond to the boundary-free and 
boundary-induced parts, respectively. It can be seen that both parts are diagonal. The latter 
property is not a consequence of the problem geometry. As it has been shown in [43] , for a scalar 
field obeying Dirichlet and Neumann boundary conditions on a plate at z = 0, the off-diagonal 
component (0 | T\3 1 0) does not vanish for a general curvature coupling parameter £. Only in the 
special case with £ = 1/4 the VEV of the energy-momentum tensor is diagonal. 



4.1 Boundary- free part 

For the VEV in the boundary-free cosmic string geometry we find the expression (no summation 
over v) 

poo poo \3 

<0 1^1 °>. = 3 E / dk dX-gfiXr), (4.3) 

j JO JO 00 

with the notations 

4%) = ^9f{y) = ~[Jj{y) + 4 + ^ 

For r 7^ the renormalization of the boundary-free part is reduced to the subtraction of the 
corresponding VEV in Minkowski spacetime: 

< T ;>s, rC n = < o i T ;! o >s-< o ! T ;! o >M- (4-5) 
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In Appendix [A] we show that the renormalized boundary- free VEV is presented in the form (no 
summation over u) 



(T v ) 

\ v I s,ren 



2m 4 



7T- 



EM 



H cos((?7r/2) / dx 



cos (irl /q) (2mr sin(7rZ /<?)) 

sinh (gx) sinh(x)-Fj°^ (2mr cosh x) 



cosh(2gx) — cos(<2"7r) 



(4.6) 



where we have defined 



4°\y) 



-3/a(i/)-/i(y), 



(4.7) 



with the notation (|3.14p . For q < 2 the first term in the square brackets is absent and the 
formula (j4.6j) is reduced to the one derived in (29j . For general case of q an alternative integral 
representation is given in [38J . For a massless field the corresponding renormalized VEV for the 
energy- momentum tensor was found in [5J [10] : 



TO, 



(g 2 -l)(7g 2 + 17) 
2880vr 2 r 4 



diag(l, 1,-3,1). 



(4.8) 



Fermionic propagators for a massive field are considered in Refs. [22, 23, 29J. Note that in the 
boundary-free cosmic string spacetime one has the relation (Tq} = (^~f) sren f° r a massive 
field. The latter is a consequence of the boost invariance along the axis of the cosmic string. 

For a massive field and for points near the string, mr <C 1, the leading term in the cor- 
responding asymptotic expansion is given by (|4.8p and the renormalized VEV diverges on the 



string as r . At large distances, mr 3> 1, for q < 2 the first term in the square brackets in 



TB]) is absent and the dominant contribution in the second term comes from the region near 
the lower limit of the integration. In this case \T^) is suppressed by the factor e~ 2mr . For 
q > 2 the dominant contribution at large distances comes from the I = 1 term in (|4.6p and the 
suppression factor is e - 2mrsm ( w /i) _ Note that the contribution of the first term in the square 
brackets of (|4.6p is always negative for v = 0,1, 3, and positive for u = 2. The contribution from 
the second term has opposite sign for 3 + 4n<g<5-|-4n with n = 0,1,2,.... However, for 
large q this contribution is small. 

In figure [2] we plotted the boundary-free parts in the energy density (full curves) and the 
azimuthal stress (dashed curves) as functions of the distance from the string in units of the 
Compton wavelength. The numbers near the curves correspond to the values of the parameter 
q. Note that the quantity (no summation over v) — (T^) STen , with v = 1,2,3, is the vacuum 
effective pressure along the corresponding direction. Hence, in the boundary-free cosmic string 
spacetime the vacuum azimuthal pressure is negative, whereas pressures along the radial and 
axial directions are positive. 



4.2 Boundary-induced part 

For the boundary-induced part of the vacuum energy-momentum tensor we have (no summation 
over v = 0, 1, 2) 



n roo roc \3 

(TZ) h = -^m(m + d 2z )J2j o d\ J dk-g^(Xr) 



cos(2kz) 
k 2 + m 2 



(4.9) 
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0.0 0.5 1.0 1.5 2.0 

mr 



Figure 2: Boundary-free parts in the energy density (full curves) and the azimuthal stress 
(dashed curves) as functions of the distance from the string for separate values of the parameter 
q (numbers near the curves). 



with the functions g^{y) defined in (14. 4p . and the corresponding axial stress vanishes: 



0. 



(4.10) 



For a massless field the boundary-induced part in the VEV of the energy-momentum tensor 
vanishes. Note that the latter is not the case for a massless scalar field (see [13]). It can be 
explicitly checked that the VEV (|4.9p obeys the covariant conservation equation V^(T^)b = 0, 
which for the problem under consideration is reduced to the single equation 



d r {r(Tl) 



We have also the trace relation 



(4.11) 



(4.12) 



Further transformations for the components of the boundary-induced part in the VEV of 
the energy-momentum tensor are given in Appendix [A] They are presented in the form (no 
summation over v) 



where 



2m 4 
^2" 



( — 1) cos (vrZ/g) F u (2mz, 2mr sm(-Kl/q)) 
sinh (qx) sinh (x 



1=0 



2q 

-\ COS (Q7r/2) 

7T 



Fo(x,y) 
F2{x,y) 
F 3 (x,y) 



dx 



cosh(2(7x) — cos(g7r) 



F u (2mz, 2mr coshx) 



Fi{x,y) = -G 2 (x,y), 
2G 2 {x,y) + F(x,y), 
0. 



In ()4.13p . the function G u (x,y) is defined as 

/•oo 

G u (x,y) = e x / due- u f v (y^T^) 

J X 



(4.13) 



(4.14) 



(4.15) 
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and the function F(x, y) is given by (|3.23|) . As we see, the boundary-induced part in the energy 
density is equal to the radial stress: (Tq\ = (Ti) h and the axial stress vanishes. By using 
the relation (|A.7|) we can see that d y [yFi(x,y)] = F2(x,y). With this relation the covariant 
conservation equation (|4.1ip is explicitly checked. 

Alternative expressions for the boundary-induced parts in the energy density and in the 
azimuthal stress are given by formulas (|A.4|) and (|A.13P in Appendix [Aj In particular, from 
(IA.4P it follows that the boundary-induced part in the energy density is always negative for 
r > 0. For the investigation of the boundary-induced VEV near the cosmic string it is more 
convenient to use the representations (IA.4D and (|A,13j) . The leading contributions come from 
the terms with j = 1/2 in the series for I(q,y). To the leading order one gets 

gm g+3 r g-i f 1)/2 (2mz) - 2mzf {q+3)/2 {2mz) 

< T A - 2[q - w (g + 2)r(( g + i)/2) ^(M^.o), ( 4 - 16 ) 

and the boundary-induced part vanishes on the string for q > 1 and z^O. Note that from the 
relation (|A.8|) in Appendix[A]it follows that the function in the numerator of ()4.16j) is negative. 

For q = 1 there is only contribution from the term with I = and from (I4.13P we obtain 
the result for a plate in Minkowski spacetime (note that in this case the asymptotic expression 
QHniD is exact): 

= ^2^ 2mz ) - 2mzf 2 (2mz)} diag(l, 1, 1, 0). (4.17) 

In this case the vacuum stresses along the directions parallel to the plate are equal to the the 
energy density. This result is a consequence of the Lorentz invariance of the problem along these 
directions. Note that in the presence of the cosmic string this invariance is broken. However, 
the radial stress remains equal to the energy density. Formula (|4.17p is a special case of a more 
general result given in |52j for a plate in flat spacetime with toroidally compactified spatial 
dimensions. 

Extracting the term I = in (I4.13p . the boundary-induced part in the VEV of the energy- 
momentum tensor can be written in the form 

( t x = + to? . ( 4 - 18 ) 

where (Tfy^ comes from the nontrivial topology of the cosmic string spacetime. We can see 
that the latter is finite everywhere except at the point r = z = 0. The divergences on the 
boundary are contained in the Minkowskian part. Near the boundary, for z < r, to the leading 
order one has (no summation over v = 0, 1, 2): 

TOb - <^>b M) - "Y^pr (4-19) 

For points away from the boundary the renormalized VEV of the energy-momentum tensor 
is presented as 

where the boundary- free and boundary- induced parts are given by (|4.6p and (|4.13p , respectively. 
The dependence of the corresponding energy density (left plot) and the azimuthal stress (right 
plot) on the distances from the string and from the boundary is presented in figure [3] for the 
geometry of cosmic string with q = 3. In the case of the energy density both the boundary-free 
and boundary-induced parts are negative. For the azimuthal stress, the boundary-free part is 
positive whereas the boundary-induced part is negative. The latter dominates for points near 
the boundary. 
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0.5 



1.5 



Figure 3: VEV of the the energy density, m 4 (TQ) ren (left plot), and the azimuthal stress, 
m~ 4 (T2 )rcn (right plot), as functions of the distances from the string and from the boundary for 
q = 3. 

5 Conclusion 

In this paper we have considered the interplay of topological and boundary-induced quantum 
effects in the fermionic vacuum. The nontrivial topology of the background spacetime is due 
to the presence of an idealized cosmic string and the fermionic field obeys MIT bag boundary 
condition on a plane perpendicular to the string. As physical characteristics of the vacuum state, 
we have investigated the FC and the VEV of the energy-momentum tensor. For the evaluation 
of these VEVs a complete set of mode functions for the corresponding boundary value problem 
is necessary. A crucial point in this paper was the obtainment of both positive- and negative- 
energy wave functions in section [2j These functions are given by the expressions (|2.29p and 
(I2.43|) . respectively. They can be used for the investigation of various quantum electrodynamics 
processes around the cosmic string in the presence of boundary involving electrons and positrons. 

The FC is decomposed into the boundary-free and boundary-induced contributions. The 
renormalized boundary-free part is given by (|3. 13j) and, for a general value of q, it vanishes for a 
massless field. As to the boundary-induced part, it is finite for points outside of boundary and is 
given by the formula (|3.22p . An alternative representation is provided by (|3.16p . The expression 
for the boundary-induced part in FC is simplified for the case of a massless field, Eq. (13.17p . In 
this case the boundary-free part vanishes and the boundary-induced part is negative for points 
outside the string. For points near the string the boundary-induced part in FC behaves as r 9 " 1 . 
The latter is the case for a massive field also. On the boundary the FC is divergent. For points 
near the boundary the part (t/S^)[ M ^ dominates and the FC is negative. The remained part in 
(V>V)b corresponds to the contribution induced by the nontrivial topology of the string and for 
r ^ it is finite on the boundary. At large distance from the boundary and for a massless field 
the leading term in the boundary- induced FC is given by (|3.2ip . Compared with the case of 
the Minkowski bulk, the presence of the cosmic string leads to an additional suppression in the 
form of the factor (r/2z) q . For a massive field and at large distances from the boundary, the 
boundary-induced part of FC is exponentially suppressed. 

Another calculation developed in this paper is the VEV of the energy-momentum tensor. 
Similar to the case of FC we have the decomposition (|4.2p . In this way, for points away from 
the boundary, the renormalization is reduced to the one for the boundary-free part. For the 



16 



general case of the parameter q, the renormalized value for the latter is given by (|4.6p . This 
formula generalizes various special cases, previously discussed in the literature. For points near 
the string the pure string part diverges as r -4 . In the opposite limit, at distances larger than 
the Compton wavelength for a fermionic particle, (1^) ren is suppressed by the factor e~ 2mr for 
q < 2 and by the factor e -2mrsin(7r/<?) f or q > 2 . 

The boundary-induced part in the VEV of the energy-momentum tensor, (T^\, is given 
by the expression (|4,13p . This part is diagonal and vanishes for a massless field. Note that for 
a scalar field with general curvature coupling parameter the part in the VEV of the energy- 
momentum tensor induced by Dirichlet and Neumann boundaries perpendicular to the plate 
has nonzero off-diagonal component (T$\ (see [S]). The latter is zero in the special case of 
the curvature coupling parameter £ = 1/4. Another difference, compared with the fermionic 
case, is that for a massless scalar field the boundary-induced part does not vanish. For a 
massive fermionic field the boundary-induced part in the radial stress is equal to the energy 
density and the axial stress vanishes. For q > 1 and the boundary-induced part vanishes 

on the string with the asymptotic behavior given by (|4.16p . Near the boundary the VEV is 
dominated by the part corresponding to a plate in Minkowski spacetime. The latter is given by 
the expression (14.171) . The corresponding radial and azimuthal stresses are equal to the energy 
density and the latter is negative. For points near the string, the boundary-free part in the 
vacuum energy-momentum tensor dominates. The corresponding energy density is negative, 
whereas the azimuthal stress is positive. 
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A Evaluation of the energy density and azimuthal stress 

In this appendix we give some details related with the transformation of the expressions for the 
boundary-induced parts in the VEVs of the energy density and the azimuthal stress. First we 
consider the energy density. Our starting point is the expression gSD with ^ 0) (A,Ar) defined 
in (I4.4I) . As the first step we insert the representation 

w = -2- r dsd s2 e-( k2+x2+m2 > 2 . (A.l) 
V 71 " Jo 

Changing the order of integrations in (|4.9|) . the integral over A is expressed in terms of the 
modified Bessel function of the first kind. After the integration by parts in the integral over s, 
we find the expression 

(T °Y = JL m (m + d 2z ) r dyJye- y 

xZ{q, y) r dk ^4 e -(* 2 WJr'A, (A . 2) 

with I(q,y) defined in (|3.1ip . The presence of the factor (k 2 + m 2 ) -1 in the integrand of (|A.2p 
leads to some complication in the further evaluation of the VEV when compared with the case 
of the FC. 
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Here we use the relation 



{m + d x )^% = e mx 
k £ + m z 



;»oo 

/ dte- mt cos(kt). 

J X 



(A.3) 



Substituting this into (|A.2p we first integrate over k with the result: 



Ctf>b 



Q m e 2mz 



7r 2 r 4 



dte 



-2mt 



d y ye- m2r2,{2y) - y{2t2,r2+l) 



l(q,y). (A.4) 



Now, by using the formula (|3.12p , the integral over y is given in terms of the Macdonald function 
and we come to the final expression 



2m 4 



(-1)' cos (nl/q) G 2 (2mz, 2mr sin(7r//g)) 



1=0 



2q f°° sinh (qx) sinh (x) . 
H cos(g7r/2) / ax — — — — r — -G2(2mz, 2mr coshx) 



7T 



cosh(2gx) — cos(g7r) 



(A.5) 



where the function G u (x,y) is defined by (|4.15|) . By using the relations 

dxU(x) = -xf u+ x(x), 
xd x f v (x) = -f v -i(x)-2vf v (x), 

for this function G u (x,y) we get the following formula 



(A.6) 



l G v+x (x,y) = (2u - 1) G v (x, y) + U-i(V* 2 + V 2 ) ~ xf»(Vx 2 + y 2 ). (A.7) 



In particular, G u (x, 0) is expressed in terms of the Macdonald function (see also [51] 

xf v {x) - fy-i{x) 



G u (x,0) 



2v-\ 



(A.t 



(2) 

Now we turn to the azimuthal stress. The corresponding expression frorngj]), with^ ; (A,y) 
defined in (|4.4p . can be written in the form 



22/b 



2m— (m + d 2z ) } ) j {2qj -l — rd r 



7T 

;°° „ cos(2fcz) 

X I dk 1,2 1 2 



oo \ 

dA ~ J i-i/2( Ar )- 
o w 



(A.9) 



By using the representation ()3.9p we evaluate the integral over A in terms of the modified Bessel 
function. Making use of the relations 



{qj - 1/2 - yd y )e y I qj -i/ 2 {v) = ye y [I qj -i/ 2 (y) ~ ^+1/2^)], 



and 



qj[igj-i/2(y) - i q j+i/2(y)} = (yd y -y + i/2)[i qj _ 1/2 ( y ) + i gj+1/2 (y)}, 

we arrive to the expression 

4mo 



(A.10) 
(A.ll) 



T 2>b = 1^ dsy 2 e-y(yd y -y + l/2)l(q, y ) 

x r dk ( m+ d 2z) c ^4e-^ 2 > 2 , 

In k + m 



(A.12) 
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where y = r 2 /(2s 2 ). The /^-integration is explicitly done by using the representation ()A,3p with 
the result 



(T|) b = -2^e 2 ^/ dUT™ I dyy 



xe -m^/2y-y(2t*y + l) {yQy _ y + 1/2 )I( q ,y). (A.13) 

Now we employ the formula (13.12p . After the integration over y we find the expression 
(T|) b = —3- J] (-l) / cos(7r//g)F 2 (2mz,2mrsin(vr//g)) 



z=o 



2q . . . f , smh (gi sinhx „ , , 

H cos(g7r/2j / ax ? — -P2(2mz, zmr cosh x j 

7r ,/q cosh(2<7x) — cos(g , 7r) 



(A.14) 



with the notation 



F 2 (x,y) = y 2 G 3 (x,y) - G 2 (x,y). (A.15) 

An alternative expression (|4.14p for the function F 2 (x, y) is obtained by using the relation (|A.7|) . 

Having the energy density and the azimuthal stress we can find the radial stress with the 
help of the trace relation (|4.12p . By using the expression (|3.22p for the FC, we can see that the 
boundary-induced part in the radial stress is equal to the energy density. 

Note that the transformations described above can also be done for the boundary-free parts 
in the energy density and the azimuthal stress. By using (|3.9p we find 



poo 

<0|^o°|0> s = dyye-^y-yi( q ,y), 

POD 

<°I T 2|0> S = 4l / dyye- m2r2 / 2 y-y(yd y -y + l/2)l(q,y). (A.16) 

Now substituting (13. 121) . we see that the contribution of the / = term corresponds to the 
VEV in boundary-free Minkowski spacetime. The latter is subtracted in the renormalization 
procedure and for the renormalized VEVs of the energy density and the azimuthal stress we find 
the representation (|4.6p . From the boost invariance along the axis of the cosmic string we have 
(0|r 3 3 |0) sren = <0|T °|0) sren . The radial stress is found from the trace relation (0|T^|0) s ren = 
m(V'V ; )s,ren) by using the expression (|3.13[) for the FC. 
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